SPECTRAL CONFINEMENT AND CURRENT FOR ATOMS 
STRONG MAGNETIC FIELDS 



IN 
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Abstract. We study confinement of the ground state of atoms in strong mag- 
netic fields to different subspaces related to the lowest Landau band. The 
results obtained allow us to calculate the quantum current in the entire semi- 
classical region B <C Z^. 
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1. Introduction and main results 



1.1. Introduction. 

Large atoms (e.g. iron) subject to strong magnetic fields exist in nature on the 
surface of neutron stars. A large amount of research in physics and mathematics 
has been devoted to the study of this system. Of particular importance for the 
present work are the articles by Lieb, Solovej and Yngvason )LSY94al ILSY94b) 
(see also |Yng91| ) which can be seen as the mathematical starting point of the 
investigation of the limits and approximating models considered in this paper and 
contain a large number of references to earlier work in the physics literature. 

We describe the atoms by non-relativistic quantum mechanics in the fixed- 
nucleus approximation. Remembering the spin of the electrons, the basic operator 
for an atom in a magnetic field is therefore the Pauli Hamiltonian 

j=l ^ 11/ i<j<fc<Ar I I 

Here Ha = (p + A(a;))^ + a ■ B(a;), with B — curl A, and a — (cti, cr2, 173) is the 
vector of Pauli spin matrices, 

= ( ? J ) ' ^2 = ( ^ V ) ' = ( -1 ) ■ 

The operator H{N, Z, A) acts on the electronic Hilbert space (including spin) 
H — /\jLiL'^{R^-C^), p = (pi,P2,P3) = (— iV) and a superscript {j) denotes 
that the corresponding operator acts on the j-th factor in the product. In par- 
ticular x'^^ is the coordinate of the j-th electron. We will also use the notation 
Pa = (pa,i,Pa,2,Pa,3) = (p + A). The scalar product in H is denoted by (•, •). 

We consider the case of constant magnetic field, i.e. we fix B = (0, 0, B) and 
A(a::) = ^B x x. The magnetic field will be strong, which mathematically means 
that we study the limit B — > +00. At the same time the atoms will be large, which 
informally means that N ^ Z and that Z — > 00 (the restriction = Z is slightly 
too strong, see Theorem 11.31 for the actual assumption) . The limiting behaviour 
depends on the relative size of Z and B. 

In the constant field case, i.e. with B, A as above, we write H{N, Z, B) instead 
of H{N, Z, A), and define the ground state energy of the atom by 

E{N,Z,B) =iniSY>ecH{N,Z,B). (1.2) 

We sometimes denote the quantum energy E{N,Z,B) by E'^{N,Z,B) in order 
to distinguish it from other energies appearing in the paper. It is known that 
for A^ < Z + 1, and all B > 0, the ground state energy E{N,Z,B) is a discrete 
eigenvalue below the essential spectrum of H{N, Z, B) (see |AHS81j V 

Heuristically one can get the right order of magnitude of the energy from the 
following description (taken from |LSY96j )^. 

• B < Z^/^ 

For small B, we can think of each electron as occupying a spherical region 
of space of radius a. The kinetic energy per particle is therefore of order 
a^^ . The small spheres organize to form a large sphere of radius R, in order 

^We use the intuitive notations <C~;^, in discussions, results will be more precisely stated. 
For instance, the statement B <C Z*/^ means that we consider sequences {(Bn , ■Zn)}neN such 
that BnZ„ — > (together with the standing assumption that Z 
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to minimize the electrostatic energy which becomes of order Z/ R. Using 
the volume relation w Na^ and setting the kinetic energy equal to the 
potential energy, we find 

The last expression, E ~ —Z''/'^, is the order of magnitude of the energy in 
standard, non-magnetic, Thomas-Fermi theory. Notice that the magnetic 
field did not enter in the discussion. 

In this region standard Thomas- Fermi theory correctly describes the ground 
state energy to leading order. 

• Z^/^ < B < Z^. 

For large B, the magnetic length scale B^^/^ becomes smaller than the 
radius of the electron 'sphere'. The shape of the 'electron' now 
becomes that of a cylinder, with axis parallel to the magnetic field, of 
radius B~^^^ and length L. The kinetic energy per particle is L~^, since the 
Pauli kinetic energy vanishes in the perpendicular variables. The electronic 
cylinders organize in a sphere of radius R. Proceeding as before, we find 
Z"2/3^-i/5^ Z^'^B-'^I^, E Z'^I^B'^/^. 

In this region a magnetic Thomas-Fermi theory (MTF) correctly describes 
the ground state energy of the atom to leading order. For B ^ Z"*/^ this 
MTF-theory simplifies since only the lowest Landau band has to be taken 
into account. 

• < B. 

When B is above Z^ the length of each individual cylinder becomes com- 
parable to the radius of the atom and a spherical arrangement ceases to be 
possible. The atom as such becomes cylindrical with radius R and length 
L with 

L ^ Z-^[\og{B/Z^)]-\ R ^ ./Z/B, E ^ - Z\\og{B / Z^)f . 

In |LSY94a] a density matrix functional was introduced and analysed and 
it was shown that it correctly predicts the ground state energy of the atom 
for B '3> Z^l'^ in particular it is valid for Z^ <B. 
For later notational convenience we define a function E — £{Z, B) that gives the 
magnitude of the ground state energy. 

rz7/3, B<z^/\ 

£{Z, B) := \ ^4/3 <B< 2^3, (1.3) 

\z\\og^)\ 2Z^<B. 

It was proven by Lieb, Solovej and Yngvason that, for all A > 0, there exists C > 
such that, if A = N/Z and Z >C, then 

C-^£{Z,B) < \E'^{N,Z,B)\ < C£{Z,B). (1.4) 

Actually they prove much more, in particular, they introduce a number of approxi- 
mating functionals, depending on the asymptotic regions above, and prove conver- 
gence of the ground state energy of the quantum model to that of the approximating 
functional. In the regime(s) where B <^ Z^ one possibility for the approximating 
model is a Thomas-Fermi type theory, depending on the magnetic field, which sim- 
plifies in the limit Z^/^ ^ B since only the lowest Landau band needs to be taken 
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into account. Magnetic Thomas- Fermi theory wiU be discussed in subsection 11.21 
below and in Section |2| For work on these questions see |LSY94allLSY94b] . |HS01| 
and |ES97[ IeSQQI IES04bl lES04a] for the case of non-constant magnetic field. 

It is of interest to determine how well these approximating theories reproduce 
the results of the full quantum model. Apart from the known approximation 
of the ground state energy, one can ask whether the ground state density and 
current are correctly predicted. The density p G L^{M.^) of the wavefunction 

G A^^^L^{R^,C^) is the function^ 

p{x) — N i \%p{x,X2, ■ ■ ■ ,xn)\^ dx2- ■ ■ dx]s . 

The current j of ■0 takes values in and is most conveniently expressed in the 
following weak sense. For all a e C^(R'^,R"^) we define 

/ a.jda;:=(7/^,J(a)V'), (1.5) 

where, with b = 62, 63) = curia, 

N 

J(a) = {a(a;(^') • p^' + p^^ • &{x^^^) + h{x^'^) ■ (t^^'^}. (1.6) 

The leading order behaviour of the density was already calculated in |LSY94al 
ILSY94b| and it was seen that all the approximating models give the correct asymp- 
totics in their respective regimes of validity for the energy. The calculation of the 
current is harder. In FouOlb it was proven that a sequence of approximating 
ground states ipz,B, i-e. a sequence of normalised functions with 

{il^z,B,H{N, Z, B)^z,b) - E^iN, Z, B) = o{eQ{N, Z, B)), 

does not necessarily give the correct current. However, the main result of the 
present paper is that MTF correctly predicts the leading order term of the current 
in its entire regime of validity. This is stated more precisely as Theorem II . II below. 

Before continuing the discussion of the current let us recall that the magnetisa- 
tion M of the atom is related to the current by the relation 

curl M = j . 

We also recall that the simpler question of calculating the total magnetisation has 
been answered in |Fou02| . 

1.2. Results on the current. 

Let T be the symmetry of reflection in the plane perpendicular to the magnetic 
field: 

{Tij){x^^\...,x^^'>) , (1.7) 

with X ~ {xi,X2,X3) = {xi, X2, —X3). Clearly, the symmetry T commutes with 
HiN,Z,B). 

^Remember that i/" takes values in -^C^ = , so the norm, | ■ j, in the expression for p, 
is the Euclidean norm in . 
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We now introduce the MTF-functional (see |LS Y94a[ IDj Y94b[ lEbiiOS] for further 
details on MTF-theory) . To a density p and a magnetic field B one can associate 
an energy Sf^ [p] by 

£n'[p]= f (r^J!^Mx))-Z^)dx + D{p,p), (1.8) 
where D(f,g) denotes the direct Coulomb energy, 



and the 'kinetic energy', T^'^^{p), is defined by 

rr^(p) :=sup(p«-n(«)), 

•u>0 

with 



, oo 

P^i-) i^i-'^' + ^&b^-f-")- (1-9) 



In the expression for Pf,, [•]- denotes the negative part: :— max{0, — x}. 

Notice that we can take a non-constant magnetic field in the definition of the 
MTF-functional. 

We denote by E^'^^{N, Z, B) the atomic ground state energy in magnetic Tho- 
mas-Fermi theory — with N electrons, nuclear charge Z and in the presence of the 
external magnetic field B — defined as 

E'^^^{N,Z,B):= inf £fl^[p], (1.10) 

with 

Cjv.B |pe i^(M^) 0<P, p{x)dx<N, 

D{p,p) <oo and / T^J^..{p{x))dx<oo]. (1.11) 

Let furthermore = (0,0,1) be the third standard unit vector in R'^. The 
analysis of S^"^^ shows that there exists a minimizer p^J^j^ z given B, N, Z. 
Furthermore, such a minimizer 'lives on the length scale t with 

^ Z"l/3(i ^5/^4/3^-2/5 0<B<Z^ (1.12) 

in the sense that the scaled density, 

Z-2(l + i?/Z4/3)-6/VMTF^^^(^,), 

has a weak limit as Z ^ oo, BjZ^ with N jZ fixed and BjZ^I'^ tending to a 
limit /3oo € [0, +cxd]. This length scale is in agreement with the heuristic calculations 
in Section [l . II fdenotcs i? there). The weak limit thus obtained is the minimizer of 
a natural limiting functional, but we will not use that fact here. 

Theorem 1.1. 

Let a^c e C^(M3,R3) ^^^^ ^^.ji^e 

a{x):^tasc{x/t), with t Z-^/^{l + B / Z'^'^y^''^ . 
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Let A > 0, /3oo G [0, +00], let (iV, Z, B) =^ {Nn, Zn, Bn) be a sequence such that 

X = N/Z, Z-^00, BZ-^'^ [3oo, BZ-^^0, 

and let ijj = ipN.z.B be a sequence of ground states of H{N, Z, B) satisfying 

Tip = ■0, or Tip = —ip, (1-13) 

then 

d 

111, .lia\ih\ — 

S{Z,B) 

Informally stated, when N k, Z and B Z^, the quantum mechanical current, 

(V,J(a)V), 

coincides to leading order with the current in MTF-theory, 

-(E^^^(N,Z,Be,+tBcmla]) . 
dt\ J t=o 



(V-, J(a)V') - 4 f E^'^^iN, Z, Se, + iScurla)) 

dt\ J t=a 



0. (1.14) 



Remark 1.2. 

The MTF-current 



-(E^^''(N,Z,Be,+tBcmlaL), , 
dt\ ) t=o 



is generally of the same order of magnitude as £{Z,B). Therefore the two terms 
on the left in (|1.14() are generally of higher order than their difference. 

Partial results on the current of large atoms have been obtained in previous work. 
It was proved in [FouOla] that MTF-theory correctly gives the current in the regime 
of non-dominant fields BZ~^/^ < C. That result was later extended in |Fou03| to 
allow for magnetic fields strong enough to confine to the lowest Landau band — 
the precise restriction imposed on the magnetic field strength being for technical 
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reasons B <^ Z st . It is the objective of the present paper to extend the validity 
of this last result to the entire MTF-region, B <C Z^ , thus proving that, when 
MTF-theory correctly predicts the ground state energy, it also correctly predicts 
the ground state current. The main improvements of the present paper over |Fou08| 
is a much more precise estimate on confinement to the lowest Landau band (given 
in Theorem 11.31 belowl . and the combination of that estimate — in the calculation of 
the current — with even better confinement estimates to slightly larger subspaces. 
Theorem |0 

1.3. Confinement. 

The word confinement in the title of this paper refers to confinement to the lowest 
Landau band. We proceed by properly defining this notion, which has already been 
informally invoked in the previous discussion. 

The kinetic energy operator in the coordinates perpendicular to the magnetic 
field, is defined by 

K^p\,i+ p\.2 + (T • B = pi 1 + pi 2 + Bas. (1.15) 

By an explicit calculation one sees that K is unitarily equivalent to a harmonic 
oscillator and that the spectrum of K is 2i3(NU {0}). The lowest Landau band (for 
one electron) is defined as the kernel of the operator K (acting on L^(R'^, C^)). The 
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projection Hq on the lowest Landau band for one electron has the integral kernel 
(see |FUFY92l) 

where x± — {xi,X2), and where ^ ^ q l ^ ' projection to the spin- 

down subspace. It is clear from the explicit expression for Ilo{x,y) (or even from 
the definition of K) that the magnetic length scale is of order B~^^'^ . This will be 
very important in the later calculations. 

We also define n> = 1— IIo. The projections Ho and n> depend on the parameter 
B (the strength of the magnetic field). We will sometimes (after scaling) need to use 
IIo, II> for other values of the parameter than the B appearing in the Hamiltonian. 
In that case we include explicitly the dependence of Ho on B in the notation by 
writing Ilo{B). 

We define as the projection in a|Li^^(K^, C^) to the space where all electrons 
are in the lowest Landau band, i.e. 

no^:=<=rf. (1.17) 
Now we can define the ground state energy for electrons in the lowest Landau band 

e2,„,«, B, . M Sp» KHiN, Z. B,n? .^jr,^^^^ ^^^^ 

Our notion of confinement is that the ground state energy for electrons restricted 
to the lowest Landau band, £^^„f, and the unrestricted ground state energy, E'^, 
are equal to leading order. Here 'to leading order' holds in an asymptotic regime 
(in B, Z) to be specified in Theorem 1 1 . 31 below . Let us introduce the parameter /3, 

13 := BjZ^I^. (1.18) 

In |LriV94a| it was proved that confinement holds under the condition /3 — > cx3. 
(It also follows from that paper that if the condition /? ^ oo is not satisfied, 
then confinement cannot hold.) The result from ,LSY94a| . though sufhcient for 
their purposes, does not include a precise estimate on the remainder term. A 
first precision of that remainder was contained in |Fou03| . Here we sharpen that 
estimate. 

Theorem 1.3 (Confinement to lowest Landau band). 
Let A > 6e given. There exists Ci, C2 > such that if 



the 



whe 



N/Z^X, Z>C2, and (3 ^ BZ^^'^ > C2, 
eIAN, Z, B) > i?Q(7V, Z, B) > (iV, Z, B){1 - Cl7^l), (1.19) 



_ rmin(r^A«+r^/35z-2A,r^/5), B<2Z^ 
'^'-\mm{B-'/^^), B>2Z\ ^ ^ 



The estimate (|1.19|l on the total energy imphes a strong estimate on the per- 
pendicular kinetic energy in the ground state. Given the kinetic energy operator 
in the coordinates perpendicular to the magnetic field, K, defined in (|1.15() above, 
we define the corresponding total perpendicular kinetic energy operator as 

N 

Under the same assumptions as in Theorem ll.3l we can now estimate the expecta- 
tion of K'^ in the ground state by the error term TZi. 

Corollary 1.4. 

Let the assumptions and notations be as in Theorem \l.Sl Suppose that -ip — 4'n,z,b 
is a ground state for H{N, Z, B). Then the perpendicular kinetic energy satisfies 
the estimate 

0< {ij, k^i;) <Cini£{Z,B), (1.21) 
where TZi denotes the error term from (|1.20|) . 
1.4. Organisation of the paper. 

The main part of the paper is devoted to proofs of confinement to the lowest Lan- 
dau band and to some slightly larger spaces. This analysis contains the principal 
new ideas. Sections contain a discussion of the current and the reduction of 
the proof of Theorem 11.11 to such confinement estimates which is our motivation 
for the present work. We recall the main steps in the calculation of the current in 
order to make the paper reasonably self-contained. If one is willing to accept the 
results of CoroUarv 1 1 . 41 and Corollarv l4.2l and mainly interested in the current, it 
is possible to jump directly to Sectional 

In Section |21 we recall estimates needed in the further analysis, mainly of Lieb- 
Thirring type. The important Section |31 contains the proof of the result on confine- 
ment announced in Theorem ll.3l In Section 0] we use the same type of analysis as 
in Section 131 to prove that much better estimates on localisation can be obtained if 
one replaces the lowest Landau band by a slightly larger space. 

Thanks. 

The author would like to thank B. Helffer and J. P. Solovej for discussions on this 
subject and comments on preliminary versions of the article. 

2. Useful inequalities 

In this section we recall the inequalities of Lieb-Thirring type that will be used in 
the proof of Theorem 11.31 Furthermore, we state some basic estimates on parts of 
the energy (kinetic or potential) that will be used as a priori input in calculations. 
We also give the Lieb-Oxford inequality which will be used in Sections [7| and |H1 

For a self-adjoint operator A such that Spec(A) n (— cx),0) is discrete, we write 
{ej(A)}j^j, with K S NU{-|-oo}, for the non-decreasing sequence of negative eigen- 
values (counted with multiplicities). The operators considered in Proposition 12.11 
below will all satisfy this assumption. 
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Proposition 2.1 (Lieb-Thirring inequalities). 

(i) Let V e L^+'^/^{W^), then 

^e,(p2 + ^) > -Ca J [V(x)t+''^' dx, (2.1) 
j 

where Cd is some positive constant depending only on the dimension d. 

(ii) The inequality H2.1() remains true with the same constant if we replace p by 

PA- 

(iii) IfV e £3/2(^3)^ ^^jg^ 

^ej(no(i/A + V)no) >-^bJ [Vix)]Tdx. (2.2) 
j 

(iv) //y GL3/2(R3)p^5/2(]g3)^ 

^e,(i/A + F)>-ABy'[y(^)]3/2rf^_^ llV{x)f'dx. (2.3) 

Proof. The estimate in was first proved in |LT76| . The extension in ^ (based 
on the diamagnetic inequaUty) can for instance be found in |Sim79| . Finally, the 
truly magnetic estimates (jml), (£3, these are |LSY94bl Theorem 2.7] and jLSY94bl 
Theorem 2.1], respectively. □ 

Remark 2.2. 

There exist generalisations to non-constant magnetic fields. See jES97l IES99I 
IES04bl IES04aj and references therein. 

Next we give, without proof, a result on magnitudes of different parts of the 
energy. This result follows from |LS Y94a[ ITSYM^ . 

Proposition 2.3. 

Let A > and be given. Then there exists a constant C > such that for all N, Z, B 
with 

N/Z = A, Z>C, 
and all normalised ground states ^ of H{N, Z, B) we have 

N 

{i^,Y.^p'i^)^^)<C£{Z,B). 
i=i 

Furthermore, if 4> & Co(K3), and B < 2Z^ then, with £ = Z-^/^{l + B/Z4/3)-2/5^ 
we can choose C such that also 

N 

\{,p,Y,zr'<l>{x^'^/m)\<c£{z,B). 

Finally, this last estimate remains true when restricted to the lowest Landau band 

N 

|(V',X!^^^'n^^V(a;(^V^)n^'V)| < C£{Z,B). 

For very large B, i.e. B ^ Z^, the Lieb-Thirring inequalities are too expensive. 
Then we need the following bound |LSY94bl 2.5 Theorem]. 
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Theorem 2.4. 

There exists a constant C > 0, such that for all Z, B > 0, 

N 

E e,(no{ifA - ^}no) > -CNz'{[\og{^ + i)r + 1). 

i=i 

For ij} S L^(R'^^) denote by G L^(M'^) the corresponding density 
^ r 

With this notation the foUowing correlation inequahty holds |L081 |. 
Theorem 2.5 (Licb-Oxford inequality). 

There exists a constant Clo > satisfying that if (z L^(M^^) is normalised such 
that IIV'IIl^ = 1, then 



■^^^ dx. 



where D{f,g) = //r3xr3 f{x)\x - y\ ^g{y) dxdy. 

3. Proof of Theorem 11.31 

The proof will be based on the Lieb-Thirring inequalities from Section |21 and on 
estimates on the commutator between Ho and the Coulomb potential. 

3.1. Estimates on commutators. 

One basic input to the proof of Theorem 1 1.31 is the following easy lemma. 
Lemma 3.1. 

There exists a constant c > such that for all f G C(M'^) with V/ € L°°(M'^), we 
have 

||[no(i),/]|| <c||v/||oo. 

Remember that the notation no(l) denotes the projection from p.l6|l with the 
parameter B = 1. Putting back the B's (i.e. performing a scaling) we get informally 
the relation 



[^o,T^]^B-'/'j^, (3.1) 



suggesting that Hq essentially commutes with the Coulomb potential away from a 
very small neighbourhood of the origin. This, as also used in |Fou03| . is our main 
new idea compared to |LSY94aj . In |Fou03| a fixed scale was introduced to bound 
the commutator by a constant, whereas in the present paper we will rather bound 
the commutator with a potential, essentially like the right hand side of (|3.1|l (see 
Lemma 13.31 for the precise statement). This turns out to give a much improved 
confinement estimate. 

Proof of Lemma WH 

The operator K = [no(l), /] has integral kernel 

K{x,y) ^no{l){x,y){f{y) - fix)). 
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We estimate the norm of K by 'Schm's Lemma' 

llifll < maxjsup / \K{x,y)\dy, sup \K{x,y)\dx}, (3.2) 

using that \f{y) - f{x)\ < \x - y\\\Vf\\^. □ 

To state the next lemma we introduce a bit of notation. First of all we will 
denote the positive Coulomb potential by V: 

V{x) ^. (3.3) 

Let /i, /2 be a smooth partition of unity on M: 

f!it) + flit) - 1, Mt) = 1 for aU t < 1/2, 

/i is everywhere non-increasing, /i(^) = for alH > 1. 

Define the cut-off Coulomb potentials by 



Define furthermore, with (x) = vT^ 



V<{x) B'/'v<{B'^^x) = -^^^^ , '"^^ (3.4) 



y>(a;) := B^/^VyiB^^'x) = ^'^^'^''"^^ (3.5) 

\x\ 



and 



Proposition 3.2. 

Let s E [0, 2]. M^it/i f/ie constant c from Lemma \S.S\ we have for all tp E L^(M'') and 
for all e > 0, 

1(^^,00 [Ho, l/>]n>^) I < icBi/2{e||C/^noVf + e-'||f/^^n>^f }. (3.7) 

We will generally apply Proposition 13.21 with the choice s = 5/4. This pa- 
rameter value assures that Ul" E L^/'^(R^) and C/^"^" E L^/2(R3), which will be 
needed for the Lieb-Thirring estimates. Recall from Proposition 12. II l lm|l that the 
Lieb-Thirring inequality in the lowest Landau band involves an L'^/^-norm of the 
potential, whereas the standard Lieb-Thirring inequality in three dimensions (|2.1(l 
(which will be applied in the higher Landau bands) contains an L^/^-norm, thus 
demanding less decay of the potential at infinity. 

Proposition l3 . 2l is an immediate consequence of the following lemma (application 
of Lemma [3.31 to the case 4>i = HqiP, 02 = n>?/;). 

Lemma 3.3. 

Let Hq be the projection on the lowest Landau band as in and let s E [0,2]. 

Then there exists a constant c (independent of B) such that for all 4>i,4>2 G L'^{M.^), 

1(01, [Ho, ^>]02)| <cSi/2||[/^ 0^11 ||[/2-.^2||. (3.8) 
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Proof of Lemma \3.!A 

By scaling it suffices to prove the lemma for _B = 1, in which case H3.8|l becomes 
[no(l),«>]</)2)| < c|| (^'^ill W^^sH- (3-9) 

Therefore, it suffices to prove a bound on (a;)^[no(l), w>](a;)^~* in operator norm, 
which we do using (|3.2(l . Writing K(x, y) for the integral kernel of the operator in 
question we find, using 

\Kix,y)\ < ±{xy{y)^-^e-^^--y-^'/'S{x, - y^) \vAy) " ^^>(^)|- 

ZTT 

We will use the simple estimates 

\Vvy{x)\<C{x)-\ 

^^<V2{rj), (3.10) 
and Taylor' formula. Thereby we find, with rj' = (rj^.O) = (y^^x^) — (xj^jXa), 
\K{x,y)\dy<C f [ {x)' {x + rj'f-'e-'^'^/^ix + trj')-^ dtdT]_L 

<2C f f\tr,'y {{1 - t)7j'f'' e-'^'^/Ut dTj^ <+oo, 

uniformly in x. The estimate with the roles of x and y inverted is similar. □ 

Armed with Proposition 13.21 we can now proceed to the proof of Theorem 11.31 
Remember that the ground state energy has the order of magnitude 

i;«ZV3/32/5 ^^2/5^9/5^ 

up to B = 2Z^, and 

for B > 2Z^. In order to keep track of this difference, we divide the proof in two. 
3.2. The MTF-regime: B < 2Z\ 

We first consider the case without electron-electron repulsion. The proof in this 
simpler case contains all the new ideas (compared to |LSY94al IFouOSj ') needed 
for the full atomic case and we consider that the choices of parameters come out 
clearer in this case. Afterwards we explain the extra argument necessary to handle 
the two-particle interaction. 

Case 1. No two-particle terms. 

In this case the Pauli Hamiltonian Hl.l|l is 

H{N,Z,B)^j:y'^'-\^\)- (3.11) 

Remember the definitions of Hp from (|l.lt)|) and that n> = 1 — Hq . For all subsets 
a C {1, . . . , TV}, we write a — {1, . . . , N} \ a and define 
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Then, with the constant c from Lemma 13.31 and all e > 0, we will prove below the 
inequality (as operators on ®j^]^L-^(M^, C^)): 

H{N,Z,B) >J2'^" (h'^ + (3.12) 



where, with the notation for the potentials introduced in (|3.3|l . (|3.4|) . H3.6|l and 
and 

and C > 1 is a constant that will be specified below. 

The inequality ^^T^ is an equality for the kinetic energies since these 

commute with the projections 11". For the potential terms the inequality follows 
from the decomposition 

11 1 

- {no(y<(x) + i/>(x))n> +n>(y<(x) + y>(x))no}, 

the inequalities (Cauchy-Schwarz) 

noy<(a;)n> + n>y<(a;)no < enoV<(x)no + e-in>F<(x)n>, 

and (c.f. Proposition 13. 2|) 

noy>n> + n>y>no = nopo, y>]n> + n>[y>, no]no 

Thus we have established (|3.12|l . 

We now estimate n"i7"n" and n^H^n" independently. Notice that if V e 
Aj^]^L^(M'^, C^), then n^ip is separately antisymmetric in the variables in a and 
those in a. 

To estimate H^H^H" we write for 77 > 0, 

77" = (1 - 27/) - ZV{x^^'^)) + f]H^ + r^H^, (3.13) 



jea 



ith 



= ^i/^^ - v^'eZB^/^Wiix^'^), (3.14) 



jea 



ith 



W,{x):^w,{B'/'x), w^{x):^i^+Cj^. 



(3.15) 



The operator Hi is an atomic Schrodinger operator (without the terms corre- 
sponding to the electronic repulsion) with nuclear charge 2Z and |a| electrons. By 
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the easy part of the HVZ-theorem we get a lower energy by adding electrons, so we 
may, for a lower bound, assume that \a\ = N. But for \a\ = N we have 

inf Spec = E'^{N,2Z,B). 

Notice that we here implicitly used the remark above on the symmetry prop- 
erties of 11" since we considered H" as an operator on the antisymmetric space 
AjgQL^(M'^, C^). We will make similar estimates without repeating this remark. 

We therefore get, since £{2Z, B) /£{Z, B) is bounded, 

inf Spec ffi" > -C£(Z,S), (3.16) 

for all a. 

For H2 we can estimate, using Proposition 12. II ( |m|l (i.e. the Lieb-Thirring in- 
equality in the lowest Landau band) and observing that wi € L'^^^' 



> -cBi/4,y-3/2e3/2z3/2 / wl^^x) dx = -c' rj'^/\^/^ Z'/^ (3^/^ Z'/' 
= -Cr,-^/\^/^p-^^^^Z-'/^£iZ,B). (3.17) 
So our total estimate on n"/f"n" becomes 

n"H"n" > n" ((i - 2tj)e^^^^^{n, z, b) 

To estimate i/" we start by observing that n>i/An> > ^n>(p^ -I- i3)n>. 
Therefore, we have to estimate 

n"ij"n" > in" (^h^ + + n", (3.19) 

where 



H?=Y.\^^Af + ^-ZV(x^'^), 

- E ^(PaV + f - e-'ZB-'/'W,{x(^^), (3.20) 



2' 



with 



(3.21) 



W,{x):^w,{B^f'x), ^,(.):=/lM+c^. 
Applying the magnetic Lieb-Thirring inequality. Proposition 12. II we get 

fff >-C f \b-^] dx = -CB-^'^Z^ t (1 - p-^f'^p" dp 

= _C"/3-9/iO£(Z,B), (3.22) 
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and, using W2 € L^^'^{ 

H2 >~C I [B- e-^ZB^'^W2{x)f_!^ dx 



> -c 



' / [e-'ZB'^^W2ix)f^dx 

= _C"g-5/2^5/2^-l/4 ^ _(^/g-5/2^13/6^-l/4 

= -C'e-'^^^Z-^/^p-'^^/'^°£{Z,B). (3.23) 
Using (pTTi^ . and (pH^ . we get 

n"ff"n" > n" (^-c[/3-9/i° + e-''^^z-^/^r'^/^°]s{z, b) + . (3.24) 

First choice of parameters 

For B < Z^^l^ we cannot profit from the 'positive term -^i? in H3.24|l . Combining 
with llTTH|) and we find 

eQ{N, Z, S) > (1 - 2,7)i?^«„„f (iV, Z, B) 



^(1+^-3/2^3/2^-3/20^-1/2) 

+ (^3-9/10 + ,-5/2^-1/6^-13/20)1^(^^5)^ (335) 



To get an optimal choice of parameters, 77, e, in 13.25|l we set the leading error 
estimates to be equal: 

^-3/2^3/2^-3/20^-1/2 ^ ^-5/2^-1/6^-13/20 ^ ^_ (3 26) 

This implies the choice 

^^^-9/35^-2/7^ e = /3-ii/70^i/2i^ (3 27) 

With this choice we get the estimate 
with 

R = O(/3-9/10) + 0(/3-9/35Z-2/7). 

This is the first error term in Theorcm ll.3l for B < 2Z^ and is the better of the two 
for B < 

Second choice of parameters 

For B > Z^"^/^ the positive term in H3.24|l dominates the negative terms in that 
equation — even for \a\ = 1. Therefore we can make a somewhat more natural choice 
of parameters, leading to the relative error term in 0(/3~'^/^) and thus better for 
large B, as follows. 
We choose 

MZB-^/^ = MZ^/^p-^/^, 7^^/3-3/5, (3.28) 

where M is a (large) constant to be chosen below. We start by analysing a bit 
differently from above. 
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The (bottom of the spectrum of the) operator is estimated as before by 
0(/3-^/i"£:(Z, B)). For H^, we apply the value of e from iHT^ and get 

m = E ^(Pa + f - M-iVF2(x(^)). (3.29) 

By scaling, -^H!^ is unitarily equivalent to the {Z, B) independent operator 

i;.^{^,<,.«,.,^_M-(/|f).C^)}, (3.30, 

where q — (— iV + B^^A) (notice that V x {B^^A) — (0, 0, 1)). Upon apphcation 
of the Lieb-Thirring inequality. Proposition 12 . II . we see that 

> ~CiM), 

where C(M) ^ as M ^ +00. Therefore, remembering the term in H3.19|l . 
for some sufficiently large (but fixed) M we find the estimate, for |5;| ^ 0, 




(3.31) 



When considering iJ" remember the parameter choices from (|3.28ll . Using the 
estimate in (|3.17|l we find 

U°'H^U°' > ~CM£iZ,B)U". (3.32) 

Therefore, i?" satisfies 

n"7?"n" > n" ((i - 3P~^^^)E^^^^f - CmP^^'^£{z, s)) n". (3.33) 

So combining H3.31|l with H3.33|l (omitting the positive term proportional to B), we 
find, 

n" (ff" + ij") n" > n" (e'^^^, + ?,(3-^'''\e'^^j - CmI3-'"''£{z, b)) n". (3.34) 

This finishes the proof of Theorem 11.31 in the case without two-particle potentials 
and for parameters B < 2Z^. 

Remark 3.4. 

If we include the positive term from (|3.31|) . and use that | -E^j^f | / 1 £^'3 1 ^ 1, we get 
the following more precise version of H3.34|l 

E n" (i?" + n" > - cp-^^^siz, s)) 

a 

+ E (i?clf-^r'^'^(^,S) + ^|a|)n". (3.35) 

\a\<N 

For B > Z^^/*" we have B > f3^^/^^£{Z, B), and therefore H3.35(l implies the exis- 
tence of a constant c > such that 

(1 - U^)H{1 O > + cb) (1 - U^). (3.36) 

This is the necessary assumption for application of the Feschbach method (see for 
instance IBFS98bl IBFS98a| ) to the present problem. So for these (rather large, 
i.e. B ^ Z^^^^) fields that approach might work. As can easily be seen from the 
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estimates below, (|3.36|) remains true, under the condition B ^ Z^^^^, when the 
electron-electron terms are included in the Hamiltonian H. 

Case 2. Full atomic problem. 

We now return to the atomic operator given in Hl.l|l . The analysis starts similarly 
to the above, only we get more terms. In particular, the inequality H3.12|l becomes 

> (1)" + ^") n", (3.37) 

a 

where, with rj^k '■= ^^^"^ ~ x^^\ 

jea 

j,kGa. : j<k 



and 



E {v{r,,u)-ie-'V^{r,,u)-Ce-'B^'^ul!\r,,,) 



j,k^a : j<k 

+ E {v{r,^.)^h-'VAr,^,)^ICe-'B'/'ul^'ir,,,)] 

The proof of (|3.37() is similar to the proof of H3.12|l and will be omitted. 

The idea behind the treatment of the two-particle terms is as follows. All es- 
timates are done using Lieb-Thirring inequalities (this was also the case before). 
The two-particle terms — typically rjf^ — each come with a coefficient 1, but there 
are N ^ Z oi them, so therefore, in total, they will contribute with a term of the 
same order of magnitude as the corresponding one-particle term — typically ZrJ^ — 
coming from the interaction with the nucleus. 

Let us start by estimating As before, we can get a lower bound by replacing 
i?^' (for j e a) by 5((Pa )^ + B). Therefore, we find 

jja^ajja > 1 jja ^^o^ + + n" , (3.38) 

where 
and 

Aj)\2 



3e-'B'/' ^ W,ir,,k)-I^-'B'/' ^ W^ir.^k)- (3.40) 

j^k^a : j <k j^a,k£a 
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In the first operator, ij" the two-particle terms come in with a positive sign and 
can therefore be neglected for a lower bound. The other operator ^2 "^iH cause us 
a bit more trouble. 

Let Ua = \a\ and fia = \cy\ = N — Ua- For simplicity of notation, wc will 
renumber the electrons so that a = {1,. a = {n^ + 1, . . . , TV}. If = 
then = 0. If not, we can rewrite as follows (using an idea of Levy-Leblond 
|LL69p for tt-q > 2 (for tt-q = 1 no reformulation of (|3.40() is necessary), 

N N 



le-^B^/\n^ ~ l)W2{r,,k) - le-^S^/" E ^2(r,v)) } (3.41) 

lea 



We can estimate n"{-}n" (where {•} denotes the operator inside { } in H3.41|l ^ 
using Proposition 12. II by 



n°{.}n" > inf-C / [B-4e-i^Bi/2^2(a;)-6e-i(n„- 1)5^/2^2 ( 



X — w) 



X- Zl] 



5/2 



dx 



lea 



> -ce-5/2z-i/6/3-i3/2"£(Z,B), (3.42) 
where the inf denotes ^'ni[2={zi,....z„ ^es.^""} ^''^^{wev.^}- Thus, 



N 



So we see by comparison with ()3.23(l that the estimate is unchanged by the inclusion 
of the two-particle terms into the operator. 

The same strategy applies to all the other terms: One writes the operator in the 
manner illustrated by (|3.41() and applies the appropriate Lieb-Thirring inequality. 
We omit the details. The error estimates are invariably of the same type as those 
given in the case without two-particle potentials. Therefore, we can use the same 
choices of parameters £,77 as before and get the same final estimate. 

This finishes the proof of Theorem II. 31 in the case B < 2Z^. 



3.3. The very strong field regime: B > 2Z^ 

B \2 



Notice that in this case we have S{Z, B) = Z^(\og^)^ and 



B > cZ' 



(log^) . (3.44) 



We keep the choice of e from the second choice of parameters from the MTF-regime, 
but choose 77 = B~^/^: 

Mx/B 



Arguing as previously, we get (|3.37|) and that is positive (for Af sufficiently large 
independent of B, Z). Using H3.44|l this becomes: There exists C > such that for 
alla7^{l,...,iV}, 

n"^"n" > csn". (3.45) 
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On the jo" we proceed as before, but with the new choice of rj. The Lieb-Thirring 
inequahty corresponding to the operator iJ^ from (|3.14|) in the case without two- 
particle terms becomes 

jja^ajja > f [b^/^ z'^Wi{x)f'^ dx > -c' . (3.46) 

As before, we get the same order of magnitude when we include the two-body terms. 
The other parts of are of lower order, and we therefore get a relative error of 
order 77, i.e. 

jja^ajja > ^£^^^(1 + CB-^/^)n°'. (3.47) 

This, combined with (|3.45() . gives the first error bound in Theorem ll.3l for B > 2Z^, 
i.e. the relative error B~^^'^. 

The second error bound in Theorem 11.31 was actually proved though not 
stated explicitly in jLS Y94a| . Here the commutation, i.e. the application of Propo- 
sition is not needed, only the Cauchy-Schwarz inequality'^. We therefore get, 
instead of H3.37|l the simpler estimate 

> ^n"(H" -l-7Y")n", (3.48) 

a 

with 

^":=E^A + E (i-36)nr-,.fc), 

jGa j,k^a,j<k 

and 

j^a j^k^a.j <k 

j^a.k^a 

We choose e = for some large M as before. After scahng and discarding some 

V B 

positive terms we find 

-3M-1 E Vir,^k)-5M-' J2 (3.49) 

Using the same arguments as before we find that there exists Mq such that for 
M > Mq the operator in {•} in H3.49II is positive as an operator on (g)jgQL^(]R^; C^) 
independently of the parameters {xk}kea- Thus, for M sufficiently large, 

n"7i"n" > 0. (3.50) 



■^Inclusion of the commutation in the argument would only improve the estimate ll.2Ui by a 
logarithmic factor. 
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To estimate Ti" we rewrite it as 

= (1 - 36) [h^I^ ZV{x^=^) + 

jea j,kea,j<k 



+ ■ 

jea 

Applying Theorem 12.41 to the last term we find 

n"H"n" > (1 - 3e)E2,M,Z, B) - CeNZ\[\og{BlZ'')f + 1) (3.51) 
> (1 - 3e)i?^Q„„f (TV, Z, B) - CeNZ-\[\og{B/Z^)f + 1). (3.52) 

Remembering the choice e — we find the relative error 

This finishes the proof of Theorem 1 1.31 in the case B > 2Z^. □ 

3.4. Estimate on the kinetic energy. 

We finish this section by giving the short proof of Corollarv ll.4l 

Proof. Define Hi{N,Z,B) H{N,Z,B) - ^K^ and let E^{N,Z,B) be the 

2 2 

ground state energy of Hi (N, Z, B) . By the variational principle we have 

i , K^^) = , {H{N, Z, B) - Hi (N, Z, B)}^) 
< E^{N, Z, B) - Ef{N, Z, B) 

2 

<e2,^^{N,Z,B)-E'?{N,Z,B). (3.53) 

2 

Furthermore the proof of p.l9|) holds with only notational change for Hi {N, Z, B). 
Thus 

£;?(iV, Z, B) > e2„AN, Z, B)il - C^7^l). (3.54) 

2 

Combining and we get ^1^. □ 

4. An estimate on confinement to a larger space 

The result on confinement in Theorem 11.31 is rather precise. However, for the 
application to the calculation of the current it is just slightly too crude. As can 
be seen from the proof of Theorem 1 1.31 the main error in the confinement estimate 
comes from the region near the singularities of the potential. The region in question 
has essentially the length scale B~^/^, which is very small compared to the other 
length scales of the atom as discussed in subsection 11.11 This suggests that intro- 
ducing a localisation away from the singular region, one can hope for a more precise 
estimate. It is reasonable to make such a localisation in the variable 0:3 (parallel to 
the magnetic field), in order for the localisation to commute with the projection on 
the Landau bands. That is the rationale behind the next result. It turns out to be 
convenient — in particular when the two-particle interaction is included — to make 
the localisation in frequency- instead of position-space. 

Let no,n> be as previously defined. Let 

_ f Z-2/5B-i/5^ b < 2Z\ 
^'-\z-^[Xog-l,]-\ B>2Z\ ^^-'^ 
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be the parallel length scale and define, for S > 0, projections on low and high 
frequency spaces 

Phf l{|p3|>5-iL-i}j Pif •= l{|p3|<5-iL-i} = 1 - Phf, (4.2) 

where 1^ denotes the characteristic function of the set f2. Clearly, phf commutes 
with Ho and n>. 

Intuitively speaking, the frequencies below S^^L^^ can only probe length-scales 
above SL. We want 5L to be larger than the magnetic length scale B^^/^ but 
shorter than the length L of the 'electronic cylinder'. Notice that 

Define the orthogonal projections Pq and P> on L'^iM.^) by 

PQ = ^a+ PM^> , P>-= Pif n> . (4.3) 

One easily sees that 

[Po,P>]=0, PoP>-0, Po + P>=l. 
By analogy with the localisation to the lowest Landau band, we also define 

N 

^^o^:-n^o^'' P>^:=l-Po^ (4.4) 

and 

E,oniM{N, Z, B) := inf Spec Po^H(7V, Z, B)P^ . (4.5) 

It is clear that phf tends to the identity as (5 — > cxd. Thus, we expect to get an 
improved confinement estimate, compared to Theorem ll.3l when 5 is chosen 'large'. 
Our result below shows that, when Z^ B, we can take S slightly smaller than 
unity and get an estimate with almost a factor of -j= of improvement over the TZi 
in Theorem ll.3l 

Theorem 4.1. 

Let X > be given. Then for all /i G (0, 1/2) there exist Ci, C2 > such that if 

N/Z = \, Z>C2, C2Z^<B, (4.6) 
and ( with L defined by (|4.1|l ) 

R-1/2 

C2^— <5<C^\ (4.7) 

then 

EconfMiN, Z, B) > E^iN, Z, B) > E^onf^N, Z, B){1 - Cl7^2), (4.8) 

where 

' \min(p-i/3,_^log^), B>2Z\ ^ ' 

Theorem 14. II implies a bound on the perpendicular kinetic energy of P^ip. 
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Corollary 4.2. 

Let the hypothesis of Theorem \4-l\ be satisfies and let ijj — z b be a normalised 



ground state for H{N, Z, B) . Then the perpendicular kinetic energy of ip satisfies 



the estimate 

0<{i>, K^P^^) < CiTZ2£iZ, B), (4.10) 
where TZ2 is defined by (|4.9|) . 
Proof of Corollary \4.^ 

Define Hi {N, Z, B) H{N, Z, B) - P^ and proceed as in tlie proof of Corol- 
lary O ' □ 



Proof of Theorem \4-.l\ 

The first inequality in (|4.8|) is an easy consequence of the variational principle, so 
we only need to prove the second. 

We define, for a C {1, . . . , TV}, 5 = {1, . . . , TV} \ a and 

:= n ^"0'^ n ^>'^ (4.11) 

We get the identities 

^ P" = 1, and ^ for a 7^ a'. (4.12) 

aC{l,...,JV} 

Proceeding as in the proof of (|3.37|l , we get the following operator inequality (where 
we write rjk instead of x^^^ — x^''^): 



H >J2p°'{Q'^ + Q'^jP°', (4.13) 

a 

with 

jea 

+ E (^^(-.^) - ^^4^V(r,.)n[^) - ^.pVX^"Vr(r,,)<))}, (4.14) 



and 



Q" = J2 {^A - ZV{x(J'>) - e-^ZV{x(^^'>) - e-^B^/^ZW2{x^^'>) 



E - r-'iB'^'W2{r,k) + V{r,,))}. (4.15) 



Notice right away that, since we will choose e <C 1 and |VF2(x)| < ^^/^^^^ (and 
C > 1), satisfies 

> ^ {^A^ - 4Ce-izy(x(^)) - 6Ce-i ^ ^(r,fe)}. (4.16) 

jea k^j 

Bound on Q. 

We will prove that for all /i G (0, 1/2), there exists cq > such that if 6, e satisfy 

^ [SLB^'^Y < Co, (4.17) 



eSBL 
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then 



pagapa > IbP". (4.18) 



We rewrite (|4.16|l as 



It thus suffices to prove that if (|4.17|) is satisfied then, for all z G M"^, 

P>^Wa - 4Ce-iZ^(x(-') - z)]pi''^ > ^Pi'K (4.19) 
12 Jo 

But (|4.19() is exactly the result of Lemma FA. 41 below. 

To assure the condition H4.17|l . we choose*^ 

Z 



e 



-{SLB^/'^YM. (4.20) 



5BV 

for some large constant M . 
Bound on Q. 

Throwing away the positive terms rj^ with j, A: € a, we can estimate, for some 
77 > 0, Q as 



>(l-377)(5](ffj^'^-Zl/(x(^-)))+ ^ 



-1 

jk 

j^a j,k(^a,j<k 



+ v{Qi + Q2 + Q3), (4.21) 

with 

Q^-H {^A - .^-^i?^/^^n(^Vi(.(^))n(^") - ^BiY.no^w,{r,,)ni^^}, 

jea k^j 

and 

^3 - E {^A - e,7-izn[^"V(xW)n[^^ - ^^n^V(.^.,)nlf)}. 

jea k^j 

To estimate the operators Q^ 's, we will use the same strategy as always: First we 
use the Levy-Leblond formula as in (|3.41l) to effectively have to estimate one-body 
operators, then we use a suitable Lieb-Thirring inequality — choosing the optimal 
one in each case from Proposition 12. II 

The term Q2 is estimated exactly as in Section O — the 11^^ 's surrounding the 
potential assuring that we can use the Lieb-Thirring inequality from the lowest 
Landau band. Proposition 12 . ll|IiH) — and we get 



P"Q2P" > -C77-3/2,3/2^3/2^1/4^ (4 22) 



^We can rewrite the choice of e as e = -^yj^{5LB^/'^)>^ ^ M. Therefore the facts that (sec 14. 7i 
and I4.6i 1 5 S> ^ J and <g B imply that the previously used assumption, e <g 1, is satisfied. 
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Notice that Qi is an atomic Pauli Hamiltonian for |a| electrons and nuclear 
charge 3Z. We can therefore estimate 

Qi > E^{\a\,5Z, B) > E^{N, 3Z, B) > -C£{N, Z, B). (4.23) 

Finally we consider Q^. Proceeding as for Q2 but using the standard magnetic 
Lieb-Thirring estimate, Proposition 12 . . we end up with 

P-Q,P->-C f \B-^f\x = -CB'/' [\i^p-^fn(!^\%2^p 

= -C'e^7j-'Z^B-'/\ 

Choice of parameters. B < 2Z^ . 

Recall that here 

L 

Write, for some i? > 0, 

r, = 7?/3-3/5 ^ RB-^/^Z^^^, 

and 

e^R^ with i?:= (5-1 {5LB^''^YM. (4.24) 



Then we find 

P"(0i + Q2 + Q3)i^" 



>-c[i+{^y'\{ii)\^yy^r^z^'^p<^. (4.25) 

We choose 



B 
Then 

and our final estimate becomes 

+ g")P" > F"(-Bco„f,hf - C5-^^{5LB^''^YB-^/^Z^^'^). (4.27) 
By and (|07|l we get 

E{N, Z, S) > -Bconf.hf - i^((5iBi/2)^B-i/5^i3/5_ (4 

V B 



This finishes the proof of Theorem 14. II in the case B < 2Z^. 
Choice of parameters. B > 2Z^. 

Here L = Z^^pog First we use the same approach as for the case B < 2Z^. 
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Then we get the estimate 

P"QP" > {(l-377)£;conf,hf 



C,;[l + (^)^S-V2(i,g^)-2 



-1] 

This follows from g^Sl, and 



3, and gl 

We can now choose e as in H4.20|l and 

B -.1/6 ,-11 B\ 
W B 



V = = S-H^f'^^^iSLB^^^rM. (4.29) 



Then 



(£)3/25l/4^-3/2(i„g^)-2 ^ (iog^)-2 < (iog2)-2, 
(1)3^-1/2 = ^ < 1/2. 

This leads to the first error bound in 14. 9|) for B > 2Z^. 

To get the second estimate in (|4.9ll for B > 2Z^ we estimate Q2 and Qs a bit 
differently. We take 

72 ]„„ _B_ 

77^ 6 °f (^£i?V2)PM. 

The argument applied to Q gives that Q3 > 0. We will use Theorem l2.4l to estimate 
Q2- To prepare for this we write 

+ JTl ^ - 3i?^/^(iV - l)W^i(r,,))}n«. (4.30) 

We now estimate each Wi{z) < -g^^, and apply Theorem l2.4l with = 1 to each 
of the N operators inside the {•}. Counting terms and using N k, Z , we get 

Q^>-CZ^[\og^f. 

This finishes the proof of the second error bound in H4.9() for B > 2Z^ and therefore 
of Theorem O □ 

5. The current 

5.1. Discussion. 

For reasons that will become clear later, we will generally impose the technical 
restriction on the test functions a in (|1.5|l . I|1.6|l that they be everywhere perpen- 
dicular to the magnetic field, i.e. 

a=(ai,a2,0). (5.1) 

The identity H5.2(l below is valid under this assumption and H5.1(l is also a cru- 
cial hypothesis for the validity of Theorem 15.11 The missing third component of 
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the current can be reconstructed from the remaining two using the spatial symme- 
tries of the Hamihonian and gauge invariance. So our final result on the current, 
Theorem ll.il does not suppose H5.1|l . 

In the regime where B is large compared to Z^/'^, the ground state ip is essentially, 
in the sense of Theorem 11.31 localised to the lowest Landau band, Ranll^ (with 
the notation from (|1.17|l '). The energy and the density can be correctly calculated 
to leading order by restricting to this subspace, but the current satisfies 

J(aX = 0. (5.2) 

Thus, one needs more than leading order information on -0 to calculate the current. 

One can consider the current operator J(a) as composed of two contributions, 
a spin-current B^h ■ a and a persistent current. The identity (|5.2II expresses 
that the sum of these operators cancel on the lowest Landau band. If we were to 
calculate the separate contributions of these two terms to the total current, the 
analysis would be much easier. 



5.2. Splitting the current. 

A first step towards the calculation of the current is to replace the operator J(a) 
by another operator having the same matrix element in the ground state tp and 
being easier to analyse. This was realised in |Fouflla| and the result is given in 
below. 

Theorem 5.1. 

Let a = (fli, a2, 0) G be given and define 

a = (-a2,ai,0), ao(a;) = a(a:;) - a(0). 
Define furthermore as the negative, symmetrised Jacobian matrix of a, 

(2dia2 d2a2 - dioi d-ia2 
d2a2 - diai -29201 -dsai 
9302 -dsai 

Define finally a decomposition of the Laplacian, 

92 92 . 92 



and the operators 

Jkin = 5]p^^Ma(:E(^'))p^^ + Ba ■ b(x(^)) - \^,_b:i{x^^^), (5.3) 



N 



N 



'fe)) ■ (a(a;(j)) 



(5.4) 



^ (.t(j) - .xW) • (a(a;«) - ^(xC^))) ,^ 
2^ |.0)_.(fc)|3 • (5.5) 

l<j<*:<A' 

Then, for any eigenfunction ip of H{N,Z,B), we have the identity 

(V', J(a)V') = {ip, (Jkin - JmT + Jdens)^)- (5.6) 
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Remark 5.2. 

The splitting of the (lower order) term A63 between Jdens and Jkin rnay seem a 
bit arbitrary at this point. The reason for including a part of this term in Jkin is 
that it is convenient to have the identity (|9.12() below. Furthermore, if one were 
to extend the analysis of the current to magnetic field strengths above Z^, this 
splitting seems to be the natural one, since in that case A^b^ is no longer of lower 
order. 

The identity (|5.6|l is valid for any eigenstate i/j for H{N, Z, B), in particular for 
the ground state. It comes from expressing J(a) — (Jkin — Jint + Jdens) as a 
commutator i[H{N, Z, B), O], where the operator O can be chosen as 

. N 

O=^^{ao(x«).p2^+pi().ao(x«)}. 

i=i 

The commutator vanishes in an eigenstate, 

[H{N,Z,B),0]i^) =0, 



and 15.(j|l follows (sec FouOla for details). This sketch of a proof gives an idea why 
the construction in FouQlbj of approximate eigenstates with 'wrong' current works: 
By perturbing the eigenstate tp a little bit, one can get an enormous contribution 
from the commutator {ip, [H {N , Z , B) , 0]ijj) without changing the energy (ip^Hip) 
very much. 

The restriction on a from (|5.1|l is forced upon us by the use of the formula H5.6|l . 
Only for a • B = is it possible to find an a corresponding to a. 

The right side of (|5.6|l is much easier to analyse than the left side. The operator 
Jdens = Z^^i '/'(^^^^Oj with (/)(a::) = _i Ay 63 (x), is just a sum of one-particle 

multiplication operators, and therefore 

(V", ■^DENsV') = / 4'pdx. (5.7) 

The operator Jint = ^ |a(3)-z(fc)|3 — ^ similar to the mterelectronic 

repulsion ^ \x'^^^ — x^*^-* |^^. It has been analysed and calculated to leading order in 
|Fou()lallRmI)8| . 

In the following sections, we will calculate each of the terms, ('0, JdensV'); 
{ip, JintV') and (?/;, JkinV') independently in terms of magnetic Thomas-Fermi the- 
ory. Before we do so, let us notice how those results will imply Theorem ll.il 

Proof of Theorem M.lX 

Using the symmetry (|1.13|l of -0 it suffices to consider a where 03 is odd in 2:3 : 

a3(a;i,a;2, -a;3) = -03(^1, 2:2, 0:3). 

But such an (0, 0, 03) with 03 odd in is gauge equivalent to a compactly supported 
(fli, a2, 0). Thus it suffices to prove Theorem 1 1.1 1 in the case where (|5.1|) is satisfied, 
i.e. 03 = 0. In that case we can use the formula H5.6|) to replace the operator J(a) 
by the three operators Jkin , Jint and Jdens ■ 
By Theorems 16.11 and 16 . 31 below we can write 

d ^MTF(^^ z, Be, + tScurla) = J^^^^ - J^Y + «s, (5-8) 



dt 
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where the right hand side is defined in (|6.8|l . In order to prove Theorem 1 1.1 1 it thus 
suffices to prove the three estimates 

(^,JkinV') = JKIN^ + o(5(Z,i3)), 
(V',JlNT^) = Jl^?'' + o(£(Z,i?)), 

Jdens^) = Julis + o{£{Z.. B)). (5.9) 

These three estimates are the resuhs of Theorems l9.1ll^?T1 and l7.ll below respectively. 
This reduces the proof of Theorem II. II to the proof of those three theorems. □ 



6. Current in MTF-theory 

In this section we briefly recaU resuhs on the MTF-functional. Details can be 
found in |LSY94b| and jFouOlal IFou03| (see also [LSTTI ILie81| for general theory 
of Thomas- Fermi-type models). We start by considering general magnetic fields, 
B e L^^M!^)- With £f]^ as defined in on the domain Cn,^ from ((rTT|l we 
have the following result. 

Theorem 6.1. 

• There exists a unique p^"^^ = pg"^^ G Cat.b such that 

Furthermore, there exists a critical particle number = -/Vc(Z, B) > 
such that J PbIv'z = min(A^, Nc). 

• The minimizer satisfies the Thomas-Fermi equation 

r{^\{p^^^)=[Vz{x)+p^^^.\xr+p]_, (6.1) 

for some unique ( chemical potential) p = fi{N, Z, B) such that /i(-/V, Z, B) = 
OforN>Nc. With the definition VeS := Vz{x) + p^'^^ * |a;|~i -h we can 
rewrite H6.1|l as 

p^T^^P^ld^eff]-), (6.2) 

and 

P|B|([Kff]-) - -r|B|(p"^^) (6.3) 

• Suppose that |B| > c > 0. Then, for all a G Cg^ (M^ , K'"^ ) , the map t h-> 
^MTF Z,'B + t curl a) is differentiahle at t — and defines a current 

jMTF f^y 

Ji^^^-adx:^ j^E^^^ (TV, Z,B + t curl a) | 

= /^{[K=ff]^P|'B|([K=ff]-)-^P|B|([K=ff]-)}dx, (6.4) 

where h — {bi, b2, 63) = curia. 

We now restrict ourselves to the case of a constant magnetic field B = (0, 0, B). 
Notice from H6.4|l that j^^F j_ j suffices to consider test vector fields a of 
the form a = (ai, 02, 0). 
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For large field strength the formula (|6.4|) is not convenient for comparison with 
the quantum current. In order to transform the expression we consider a = 
(ai,a2,0) and define 

a= (-a2,ai,0), ao(x) = a(a;) - a(0). 
Let furthermore, p^TF ^j^g minimizer of f and define (for small t) 

At{x) = det(/ + tDao{x)), pt{x) = At{x)p^'^^{x + tao{x)). 

Notice that At = 1 + ttr[L'ao] + ©(i^) = I - tbs + 0{t^). Define finally, the 
diffeomorphism (pt by (f)t{x + tSLo{x)) = x. Changing variables in the integrals we 
can calculate 



^t{4>t{y)) 



' Ay)dy+UfTT4^^^dxdy. (6.5) 



\My)\ ^JJ \Mx)-My)\ 

Using that the derivative of [pt] at t = has to vanish, combined with the 

Thomas- Fermi equation (|6.3() . we get the relation 

where 

1 f f ^^^{x - y) ■ (a(x) - a(y)) ^ ^ 

D-M.g) ■■=^JJ fi^y- — 1^ _\^j3 ^-^9{y)dxdy. 

Therefore we can write the MTF-current in the case of constant magnetic field as 

• MTF Tj„ ^ _ tMTF tMTF , tMTF /f. r,\ 

J • rsaaa; — jj^jN — ^int +'^denS' lO-'y 



where 



^Kl7 = / b,{[V,sUP'^mf,y)^^PBms\-)}dx, 
tMTF _ n /'^MTF „MTF\ 



Remark 6.2. 

For weak magnetic fields the discrete sum in the definition of Pb (see (|1.9|l ') can be 
approximated by the corresponding integral and one finds Pb{v) ~ const x v^^"^. 
On the other hand, if B is strong then only the first term in the sum contributes, 
and Pb{v) ~ gfj-w"^^^- For the atomic MTF-problems these approximations are 
correct to leading order if _B <C Z'^^'^ (weak field) and B ^ Z'^^^ (strong field). 
The original formula for the current, H6.4|l . thus suggests (and a rigorous analysis 
confirms) that 

J j^TF .Badx^ o{E^^^{N,Z,B)), 

for B <C Z^^^. Similarly, the changed homogeneity for large fields gives that 

J^^^^o{E^^^{N,Z^B)), 

for B > 

29 



Theorem 6.3. 

Let B = 5(0,0,1) and B,N,Z > 0. Then the MTF-current, j^'^^ defined in 
Theorem \6.1\ satisfies 

• //a II B, then /jMTF . ^ q 

• //a±B, thenj^'^^ satisfies (|^ . 

Furthermore, for all a.sc G C5"(R'^), all e > and all X > there exists a constant 
C > such that, with a(a;) = iaLsc{x/e), (i being defined in ifTT^ ) i/B < C-^Z^/^, 
ond N/Z = X, then 



J j^^^-Badx\ <e£{B,Z). 
Also, if B> CZ^/^, and N/Z = X, then 



J, 



MTF 
KIN 



< e£{B,Z). 



(6.9) 
(6.10) 



7. Calculation of Jdens 

This calculation was already carried out in |LSY94b| but for completeness we 
give a sketch of the proof. 

Theorem 7.1. 

Let slsc ~ (ai,sc, 0-2. sc, 0) G C^(]R'^) and define a{x) = £slsc{x/C) with I from H1.12|l . 
Let a. So and b be as defined previously. Furthermore, let the operator Jdens &e 
defined by 1)5.4(1 and let us define 

^■ao(x) 

^ ' |x|3 

Let A > and (3oo G [0, +oo] and let (N,Z,B) = {Nn, Zn, Bn) be a sequence with 
Z ^ OO, and such that 

X = N/Z, -^/3oo, B/Z^^O, 

(as n ^ oo). Then, if tjj — ipN.z.B *s o-n associated sequence of ground states of 
H{N,Z,B), 

^ ' ^ ' ^ ZW{x)p'^JJx)dx 



(-0, JdensV') 



0. 



(7.1) 



£{Z,B 

Sketch of proof of Theorem \7.1\ 

Since ^ Z/l for B Z'^ , the term with A||63 in Jdens is clearly of lower order 
and will not be considered. 

Let M > and consider Wm^^) = ^{\x\<Mt)^ {^)- For M sufficiently big, we 
have 

\Wm{x)-W{x)\< ^ 

and therefore. 



M2£ 



{i,,Y^Z[WM{x^'^) - W{x^^^))i,) <^lpdx< j^S{Z,B). 



(7.2) 



A similar inequality holds on the Thomas-Fermi side and it therefore suffices to 
prove 



p'^{x)WM{x)dx- Z I p^^^{x)WM{x)dx ^o{£{Z,B)), 
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(7.3) 



for all M > 0, where p'^ is the density of the ground state ip- 

Consider, for a e M, the self-adjoint operator (on the electronic Hilbert space 

n) 

N 

H^{N,Z,B) := H{N,Z,B) + aZ^ZWM(.x^^^), 

and define Ea{N,Z,B) inf Spec iJa(7V, Z, B). 

By the correspondence MTF — quantum mechanics proved in |LSY94b| . we have 
for a = 0, 

E{N,Z,B)=£fl^[p^'^^]+o{E{Z,B)) 

= -J PB{[VM-)dx - D{p^'^'' ^p^^"") - ^iN, + o{£{Z,B)). (7.4) 

We need a similar lower bound on Ea{N, Z, B). The desired bound was already 
given in LSY94b , but we recall the main line of reasoning for later reference. 

Lower bound on Ea{N, Z, B) 

For (j) e AjLii^(K^,C^) with \\(f>\\ = 1 and density p^, we estimate 
{(b, H^iN, Z, B)dp) > inf Spec l)iN,B, V,s + aWu) 

+ E iJirhm^^ ~ 2D{p^^\ p,) - f^N, (7.5) 
j<k I I 

where i){N, B, Vcs + aWM) is the mean field hamiltonian 

N 

[)(iV, B, Kff + c^Wm) = ■^-ff + ^Wm)^'\ 

i=i 

t)\B, Kff + aWM) = pi + B • + Kff + aWM- 
We now use the positivity of the Coulomb kernel 

D{fj)>0, for all /, (7.6) 
together with the Lieb- Oxford inequality, Theorem 12. 51 

(<^' E M,)!^W| ^) ^Dip„ p,) - Co I P'l' dx. (7.7) 

j<k ' ' 

So we find 

(<^, H^{N, Z, B)4>) > inf Spec \){N, B, V,s + oWm) 

- D{p^^\p^^^) -pN^ C^o j pf dx. (7.8) 

At this point we need the semiclassical asymptotics of the mean field operator 
(cf. |L^V94bl Theorem 3.1]). 

Theorem 7.2 (Magnetic semiclassics). 

Suppose that C > and that ub.z *s o, potential depending on the parameters B, Z 
and such that with h = £-1/2^-1/2^ 5 ^ ^£3/2^-1/2 . 
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• The quantity 

hb f, ,.c,/2 , 1 



l + hb 

is bounded uniformly for B < CZ^ 
For 
that 



• For all e > there exists R > independent of B, Z for B < CZ^ such 



^ [uB,zf-^ dx+ I [uB,ztl^ dx < e. 



1 + /l6 7{|3;|>fl} ' 1 + 7{|2;|>fl} 

Then, with Ub.z{x) — ZI^^ub,z{x/ I) and Pb being the pressure function intro- 
duced in p.9|l . for all rj > there exists Hq > such that if h < ho then 



tiii)\B,UB,z)]- _^ 



J^^Pb{[Ub,z]-) 



< 77. (7.9) 



Notice that when Z — > cx), then h 1 iS B Z^ . Therefore, using Theorem l7.2l 
and the Thomas-Fermi equation (|6.3|l . (|7.8|) becomes 



{^,H^{N,Z,B)cj,)>^{l + 0{l)) J PB([^eff + aV7M]-)dx-7^(p^TF^^MTF) 

-^iN-Clo J p'^J^'dx. (7.10) 

For (p such that {(f), Ha{N, Z, B)(l)) < (the only ones where a lower bound is 
non-trivial) we can estimate 



J p'/'dx<CZ-'/'£{Z,B), 



by the Lieb-Thirring inequality (see |LSY94b l p. 121] for details). The first term 
on the right hand side of ()7.10fl has order of magnitude £{Z,B), so we can rewrite 

(frrnji as 

(0, HJN, Z, 5)0) > - y" PBms + aWM]-)dx - Z?(p^TF^ ^mtF) 

- fiN + o{£{Z,B)). (7.11) 
This is our lower bound on Ea{N, Z, B). 

Finishing the proof of Theorem \7.1\ 

We now combine ()7.4|) and (|7.11|) . The terms proportional to /i cancel, since /i = 
for N > Nc. Therefore, we can estimate as follows, using ()7.4() and (|7.11() to get 
the last inequality 



aZ I p'^{x)Wm{x) dx = (V-, H^iN, Z, B)iP) - H{N, Z, B)^) 

> Ec,{N, Z, B) - E{N, Z, B) 

> J PBms]-)~PBms + aWM]-)dx + o{£{Z,B)). (7.12) 
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Applying H7.12|) for positive and negative a we find 
a > 0: 

P^{x)Wm{x) dx > Z / — —dx 



a < 0: 



-oi£{Z,B)). (7.13) 
a 

Z \ p^{x)Wm{x) dx < Z / — —dx 

+ -o{£{Z,B)). (7.14) 
a 

Remembering the Thomas- Fermi equation p^^TF _ P^[[Vcs]-), we get (17.311 from 
(fTBjl and GUI). This finishes the proof of Theorem O □ 

8. Calculation of Jint 

The calculation of {tp, JintV') was carried through in [FouOlal iFouOSj . For con- 
venience of the reader, we give an outline of a proof. 

For a function a : ^ M'^ we define the mean field interaction term Z?a as 
follows 

The comparison between (-0, JmT'4') and Jj'^x^ is contained in the next theorem. 



^a(/,.g) 77 // fi^)- \, 3 —9{y)dxdy. (8.1) 



0. 



Theorem 8.1. 

Let Bsc — (ai,sc) 12,80, 0) G C^(R'^) and define a{x) = £asc{x/i) with i from p.l2|l . 
Lei a fee as defined previously. Furthermore, the operator Jint is defined by (|5.5|) 
and by (lO) . 

Let A > and /?oo G [0, -f cxd] and let (N,Z,B) = {Nn, Zn, Bn) be a sequence with 
Z — !■ oo, and such that 

\ = N/Z, B/Z^^^^^f]^, B/Z^^Q, 

(as n ~* oo). Then, if ^ — ipN,z.B is an associated sequence of ground states of 
H{N,Z,B), 

1 

S{Z,B) 
Sketch of proof Theorem \H.l\ 

The strategy of the proof is similar to the one of Theorem 17.11 We consider, for 
a € [—ao,ao], and ao > sufficiently small, the self-adjoint operator 

K"*(iV, Z, B) H{N, Z, B) + aJm^, 

and define E'^\N, Z, B) inf Spec H'^'^N, Z, B). 

When going through the steps of the proof of Theorem 17.11 — but for the new 
operator — the main difficulty is to make sure that analogues of (|7.6(l and H7.7|l 
hold. More precisely, we need to be able to choose ao sufficiently small that 

DifJ) + aD~MJ)>0, (8.3) 
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for all \a\ < ao and all / with D{f, f) < oo. Furthermore, we need the Lieb-Oxford 
type inequality in l|8.5|l below. 

These two crucial estimates are the results of Lemma [8. 31 and Lemma [8 . 21 below. 
With these extra ingredients the proof follows that of Theorem 17.11 with mainly 
notational differences and will be omitted. □ 



Lemma 8.2. 

Let a and a^c be as defined in Theorem and let he defined by (|8.1|l . Then 
there exists a constant C > (depending only on a finite number of seminorms 
||5"asc||Loo(R3)j, such that 

\DMJ)\<CDifJ), (8.4) 
for all f e C5^(R3) with D{f, f) < oo. 

Proof. By scaling it suffices to prove (|8.4II in the case £ — I, i.e. for a = Sigc- Let JC 
be the operator with integral kernel 

^, ^ _ jx-y)- {ajx) - a.{y)) 
i^\x, y)— I |„ 

\x - y\-^ 

Recall that (47r)"^ jx — ?/|~^ is the integral kernel of the operator (— A)~^ which we 
will denote by Using integration by parts we therefore find 

47r/C = V • \asc - slsc ■ iv. 

We introduce a factor of \p\^^ on each side and get after commutation 

4-^ = • asc - a.. . V + ^ . [a,„ b|] - [b|,a,,] • 

- p{(diva.c) + ^ • [a.., H - [H,a..] . ^}^. 

To finish the proof of Lemma 18.21 we therefore only need to know that the com- 
mutator is bounded for functions (j) that are smooth and have bounded 
derivatives. This well-known fact can for instance be seen by splitting \p\ in a 
smooth, unbounded part — for which pseudodifferential calculus gives the result — 
and a compactly supported part for which the commutator is bounded as a com- 
mutator between bounded operators. □ 



We now give modified correlation inequality in the spirit of Theorem 12.51 For 
ip e L^(M^^) denote by € L^(R^) the corresponding density 

^ r 

p^{x) ■.= y2 / \i/j{x^^\ . . . ,x'^^'>)\'^S{x - X 
Lemma 8.3 (Modified Lieb-Oxford inequality). 

Let a and slsc be as defined in Theorem and let be defined by (|8.1|l . Then 
there exist constants aQ,Ci,C2 > (depending only on a finite number of semi- 
norms \\d°'asc\\L'^(R'>)), such that for all normalised ip £ L^(R'^^) we have the 
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inequality 



j<k ' ' 



> 



(l - Cia^)D{p^,,p^) + aD^{p^,,p^) ~ C2 I p^/^ dx, (8.5) 



whenever \a\ < ag. 



Proof. 

By scaling it suffices to prove H8.5|l in the case £ — 1. We choose ao sufficiently small 
that X i—> X + aa.sc{x) is invertible on M.^ for \a\ < ao. Let 0q be the inverse, i.e. 
(j)aix + asisc{x)) = X and define Aa{x) = det(l + aDsisc{x)). By Taylor's formula 
there exists C > such that 

{x - y) ■ [x - y + a{asc{x) - asc{y))] ^ o? 



\x-yf \x-y\ 
1 



< 



\4>c,{x) - 0a (y) I 



^ {x - y) ■ {x - y + oi{a,^(x) - a,c{y))\ ^ ^^^^ 

V-y\' 



Therefore, we get (for sufficiently small a), 

_a;(fe)) . [(x(j') -xW) + a(a,,(a;(j)) -a,,(a;W))] 

> (1 - Ca2)(V.,5] — — -— V') 

^ |(^a(a;W)) - ?!)a(a;W)| 

= (1 - 5: , (8-7) 

where (with the product being over 1 < < A^), 

V'o(xi, . . . , xjv) = yilM^ + aa,e(x(i)), . . . , + aa,,(x(^))) . 

Now the standard Lieb- Oxford inequality. Theorem 12.51 followed by (|8.6(l imply 
that (with Pa being the density of ^q) 

^"^"'Xl r7?r~~~(Tri^") - ^(^"'^") ~ / p'^J'^{x)dx 



j<k 

> 



> (1 - p) + ^) _ 2Clo / c^a^, (8.8) 



where the estimates are for small a. Putting together (|8.7|l and 1)8. 8|) and using 
Lemma 10 we get (|8.5|l . □ 
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9. Calculation of Jkin 

9.1. The result. 

The analysis of (■0, Jkin^) is very different in the two regimes B < Z^^^ and 
B > Z'^/^. This reflects on the one hand the localisation to the lowest Landau 
band in the high B regime, and on the other hand the vanishing of in the 

same parameter domain. The case B < Z^l'^ was already treated in |Founia| . We 
give the main ideas for completeness. 

Theorem 9.1. 

Let asc = (ai,sc, a2,sc, 0) G C^(R"^) and define a{x) — £a.sc{x/£) with £ from (|1.12|l . 
Let a, Ma, Jkin be as defined in Theorem \5.1\ and let he as defined in (|6.8|l . 

Let A > and Poo € [0, +oo], and let {N,Z,B) = [N^, Z^ B^) he a sequence 
with Z — > cx), and such that 

(as n ^ oo). Then, if ijj = 4'n,z,b is an associated sequence of ground states of 
H{N,Z,B), 

1 



(tA, JkinV') - Jkin ^0. (9.1) 



£{Z,B) 

In the case /3oo = +oo, H9.1|l is improved to 

^^^{[(V-, JkinV-)! + - 0. (9.2) 

9.2. Case of /3oo < +oo. 

Here we will prove (|9.1|l in the case B < Z^/^. We introduce 

Hl''\N, Z, B) := H{N, Z, B) + oJkin, 

and proceed as for the calculation of Jdens ■ Clearly the crucial point is to establish 
a semiclassical result similar to Theorem 17.21 This was obtained in jFouOlaj from 
which we get 

Theorem 9.2. 

Let Ub,z he a potential satisfying the hypothesis from Theorem \7.<H Define 

Pb^uAv) ■■==rJ2 d,ir^[{2,. + 1)BV, - B{1 + 2tb,{u)) - vff, (9.3) 
with do := d^, :— tt^^ for v >1. Here 



bu,t ■■= IcurUv^l + iAfa(w)A(v/l + tM4u)x) \ = 1 + tbsiu) + 0{t'^), 

Au^t := I dct{^TTtMju))\ = 1 - tb-siu) + 0{t^). 

Let A > and [3ao G [0, +oo), and let (TV, Z, B) = {Nn, Z„, Bn) he a sequence with 
Z oo, and such that 

A - N/Z, B/Z^/^ ^ /3oo, 

(as n ^ oo). Then 

-tT[l,l{B,UB.z]->- [ PB.uAUB,ziu))du + oi£{B,Z)). (9.4) 
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Using Theorem 19.21 instead of Theorem 17.21 one obtains H9.1|l for (3oo < +00 by 
the same method as for the proof of Theorem 17.11 We omit the details. Notice 
though, that the second order (in t) difference between b^^t and 1 + ^63 becomes a 
dominant term for B '3> Z"^^^. On a technical level, this is the reason why we have 
to treat the case of large B differently. 



MTF 
KIN 



9.3. Case of /3oc = +00. 

In this region we will prove the improved estimate (|9.2|l . Since, 
o{£{Z,B)) in this regime, we only have to prove the corresponding bound on 
KV'i >/kin'0)|. We state and prove a slightly more general bound. 

Theorem 9.3. 

Let hsc — {bscA,bsc,2,bsc,3) G Cq(]R'^,M'^) and define h{x) :— hsc{x/£) with £ = 
^-1/3(1 _|_ B/z4/3)-2/5_ IpJ- ^ {Msc,t,j}lj=i be a C^-function with values in 
the symmetric 3x3 matrices and such that 

A^sc,3,3 = trM,c = 26,c,3- 

Define M{x) := Msc{x/£). Let X > be given. For all e > 0, there exists C > 
such that if ijj is a normalised ground state of H(N,Z,B) where 

N/Z = X, Z>C, CZ'^I^ <B <C-^Z^, (9.5) 

then 

\{-^,JKm^)\<e£{Z,B), (9.6) 

where Jkin is the operator introduced in H5.3(l for the given M, b. 

We will for shortness write H9.6|l as (■(/', JkinV') = o{£{Z,B)) without explicitly 
including the limits (large Z, B) and dependence on parameters from Theorem l9.3l 
in the notation. We will use this shorter notation in the proof below. 

Proof of Theorem 

We will reduce the proof of Theorem 19.31 to the estimates on confinement. Theo- 
rem ^31 and Theorem 14.11 
We write AI as 

M ■.= M+JV 

hi M12 0\ / A^i3\ 
h2 M22 + A^23 . (9.7) 
D 0/ V^i3 ^23 / 

We split Jkin accordingly 

N 



2 



Ji:=PAMpA + Ba-h--Aj_b3, J2 ■= PaMpa- (9.8) 

The part of Jkin which is hardest to estimate is ^ J^^'^ . For this part we need 
very precise estimates on the confinement to the lowest Landau band. 
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Lemma 9.4. 

Let the assumptions be as in Theorem \y.y\ and let Ji , J2 he as defined in (|9.8|l . 
Then 

N 

(V',^JpV)=o(£:(Z,S)), (9.9) 

N 

(V',5]J^-''V) (9.10) 
i=i 

Clearly Theorem 19.31 follows from Lemma We prove the estimates (|9.9ll and 
(|9.1U|I separately. 

Proof of 

Reduction to the non-diagonal part. 

Notice that the assmriptions imply that Aj^^a = i~'^U{x/t) for some U E Cq(M^) 
and that for B > Z^l'^ we always have the relation < B. 

The operator K defining the Landau levels is unitarily equivalent to a harmonic 
oscillator. We can define the corresponding raising and lowering operators a* and 
a by 

a :=Pa4 - *Pa,2, a* := Pa,i + jpA,2- (9.11) 

In particular, we have alio = 0. 

A direct calculation, expressing the pA in terms of a, a*, gives 

Ho JiHo = 0. (9.12) 

Here we used, among other things, the identities 

noPAMpAHo = |notr[M]no + ino(A^tr[x])no 

Hob • aHo = -nofcaHo. 
Furthermore, there exists a constant C > such that 

±Ji < CiJi^B). 

Thus, 

±n>Jin> < 2ck. 

Using (a weak version of) Corollarv ll.4l we therefore get 

(^,5]4^')jW4^V) - o{E{Z,B)). (9.13) 

Combining H9.12|) and H9.13|l we estimate the diagonal part of ^ Jp-* . Thus only 
the off-diagonal part, 

N 

Ji^^ff := ^ j(^)4) + 4) Jp)4")), 
i=i 

remains to be estimated. The way we estimate Ji,off will depend on the magnitude 
of B. 



Case 1. < B < Z^"^!^. 
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For < B < the error bound in Theorem 11.31 is sufficient for a direct 

attack. 

The Cauchy-Schwarz inequahty and the estimate ||a*no|| < \J2B imply, for any 

7y>0, 

±(no Jin> + n> JiHo) < c^-^B + ^k. 

Taking ?/ = ^T^i^ , for some e > and with TZi defined by (|1.2U|) . we get 

KV'.Ji.ff^)! <C6-l|7^l(^)^/^{zV3(_|_)2/5} + ,7^^^(^,i^^^>. (9.14) 

Notice that, in the parameter regime studied, we have the relations 
z7/3(B/Z*/3)2/5 ^ g^z, B), 7^l « (B/Z4/3)-3/5_ 

Therefore, due to Corollarv ll.4l 19.14|l implies 

(V-, Ji.offV-) =o(f(Z,S)). (9.15) 

Combining H9.15|l with (j9.12|l and H9.13|l finishes the proof of ()9.9|l in the case 
Z4/3 « B « Z13/6. 

Case 2. Z"/^ < ^_ 

For < B we do not have a sufficiently precise estimate on the confinement, so 

we need to use also Theorem 14.11 The analysis below is valid ioi <^ B < 2Z^. 
We decompose Ji^off, as 

Jl.off = 3l +32 +53, (9.16) 

with 

Phf ■■= l{|p3|>i5-iL-i}> -Pif := 1{|P3|<|5-1L-1}: 

and (with Pq , P> and phf , Pu as in Theorem 14. 1|) 

w 




33 E (Pl/^n« Jp%[^"Vif +pLf n[^") Jp%(^")p(/)), (9.17) 

The last component, satisfies the estimate 

±^3 < CNBiSL)^, 

for all M € N. This follows by standard semiclassical pseudodifferential calculus 
as in |Rob87j . since puphf = 0, so therefore the operator in question has vanishing 
symbol. 

We estimate 3i for any ei > 0, using the Cauchy-Schwarz inequality and Corol- 
lary 

|(^,aiV'>| < c(e^'7^2A^B + el7^^^(V,i^'^P> V-)) 
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Choosing ei = {j^^Y^^ with 7^2 from gSI), we get 

/ TZ-yNB \ 1/2 

<C(^^|-^) £{Z,B). (9.19) 

Remember that 7?.2 depends on a parameter 5 which we will choose when we have 
estimated Z2- For any 62 > we get, using Corollary 1 1.41 

N 
N 

< C•(e2"^7^lS(52i2^e2)f(Z,B). (9.20) 
We choose £2 = SL^TZiB and get 



Remembering that for B < 2Z^, 



we get 



^^^^ ^ :(5Li?V2)P < ^ ^ - -5. 



£iZ,B) sVb' ' - (<5VB)i-a'' 

Comparing (|9.19() and H9.21|l . and remembering the condition (|4.7|l . we therefore 
have to choose (5^1 and /i < 1/2 subject to the restriction 

max(r3/io^^i/(i-M)5-i/2) 

Clearly this is possible in the parameter regime B > Z^+i/io^ for instance the 
choice 

^ = ^, ^ = { max (r^/^°, zV(i-M)5-i/2) 1^/^^ 
works. So we have proved that for Z^+i/i" <B<2Z^, 

{^,Ji.o«^) ^o{£{Z,B)). (9.22) 
Combining I^JE^i and (l??^ we get □ 

Remark 9.5. 

For B>2Z^. we find 



^2A^g n^B r— B 

f(^^z^(iog#.)^' '5^V^ = ^V''i^- 

We need both these terms to be 0(1). We insert 7^l = B^^/^ and 7?.2 = (5"^^S"i/3 
and find that both terms can be made o(l) as long as i? < Z'^~^^ for some /i > 0. 
Therefore, our results actually also permit an analysis of the current for magnetic 
field strengths B much stronger than Z^ , i.e. B < Z^~'^. For Z^ < B, however, 
MTF-theory fails to correctly approximate the behaviour of the ground state en- 
ergy and it is unclear what to substitute for ^ \ ^_^E^'^^ {N, Z, Be^ + tBcurla) in 
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(|1.14|l . In |LSY94a| a density matrix functional is analysed, but it is not known 
how to correctly generalise that functional to non-constant fields. 

Proof of 

Using the raising and lowering operators from H9.11|l we find (with Njk being the 
entries of the matrix M in (|9.7|l ) 

U0J2U0 = Ho [a, iVi3 + iN23]IloP3 + (Ho [a, + iiVsslnoPs)* 

We can write the commutator as [a, N13 + iN23] = l^^(t>{x/ 1) for some (f) G Co(M''). 
By the Cauchy-Schwarz inequality we therefore get 

±no JaHo < e-ir2no|0(a;/£)|2no + epl 

Since ^ Z for B <C we therefore find, using Proposition 12. 31 

N 

We also write 

n> J2n> = n>An>p3 +P3n>A*n>, 

with A being the operator 

v4:-pA,iiVi3+PA,2iV23. (9.23) 
Notice that, since the functions Nij are bounded, 

n>An>A*n> < n>^^*n> < ck. 

Therefore 

±n>j2n> < epl + €-^Ck, 
for any e > 0, which together with Proposition 12.31 and (a weak form of) Theo- 
rem ^31 implies 

N 

[i,,Y,U^i>j!i^Il'i>^)^o{E{Z,B)). 
Finally the off-diagonal terms which we write as 

noJ2n> + n>j2no = 2{n>AnoP3+P3noA*n>} + noin> + n>i*no, (9.24) 

with A from ifH?^ and 

A ■■= PA,l[^13,P3] +PA,2[^23,P3]- 

We estimate, using ||noa|| < C^/B, 

(noA*n>)*no^*n> < ck, 

UoAUy + UyA*Uo < e'^BUy + eI-^Uo\(f>{x/£)\'^Uo, 
for some function if) £ Co(R'^). Thus, applying the Cauchy-Schwarz inequality to 

±(noJ2n> + n> j2no) < epj + 2Ce-^k + er^UQ\(j){x/e)\^Uo. (9.25) 

Since Z for B <^ Z^, we therefore conclude from Proposition 12.31 and (a 

weak form of) Theorem 11.31 that 

N 
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This finishes the proof of (|9.1U|) and therefore of Lemma [9.41 □ 

Clearly Lemma 19.41 was all that remained to be proved in order to establish 
Theorem E31 □ 

Appendix A. Estimates on pu-^pu 

In this section we consider the operator 

1{|P3|<7} |^1{|P3|<7} 

for any 7 > and s > 1. We will in particular be interested in the case s = 
1. Of course, the index '3' denotes the third component; x = {xi,X2,X3) and 
p — (pi,P2,P3) ~ — «V. Informally speaking, the cut-off in frequency prohibits 
localisation in space on length scales shorter than 7"^, which explains the main 
result in Lemma lA. 31 below. 

Let / G C^(R) be an even function, satisfying that / = 1 on [—1, 1], supp/ C 
[—2, 2] and define 

f,{t) = /(7-4). (A.l) 

For reference, we first state without proof the following result of an elementary 
calculation. 

Lemma A.l. 

Let f e C5^(]R) be even. Then the operator /^(pa) as an operator on ^^(R) has 
integral kernel 

Kyixs^ys) —filix - y)), 
Zn 

where f G 5(R) is the Fourier transform of f . 
Lemma A. 2. 

Let f £ C^(M) be even. Then, for all s>l, q> 1 there exists C > such that for 
all a, 7 > we have the estimate 



1 



< c 

e(L2(R,3)) 



(a7) 



1/9 



Proof. The operator has integral kernel 

Qi{x,y)^ I K^{x,w) ^ K^{w,y)dw. 
By symmetry and Schur's Lemma, we have 

K^ix,w)\ ^ \K^{w,y)\dwdy, 

so by the Holder inequality, we get with q^^ + (g')^^ = 1) 

[ f l/(7(^-^)) l 



IQ7II < sup 



^ [ I ) I 

< 1 llfll llfll («7)'/n f du' ^1/9 
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Lemma A. 3. 

Let f G C^(R) be even. Then, for all s>l, q> 1 there exists C > such that for 
all J > and all z G R'^, we have the lower bound 

1 7^/' 

-fiiPa)-. (7/7(^3) > -C- [7^177. 

\x — z|* \x± — z±\'' '-'1 

as operators on L'^{M.^). 
Proof. 

By implementing the unitary scaling x 1— > 'y~^x and translation on M^, it suffices 
to consider the case 7=1,2 = 0. By Lemma FA. 21 

>-C f \i^{x)\''\x^\-^-' dx. 
This finishes the proof of Lemma FA. 31 □ 
Lemma A. 4. 

Let P> be as defined in (|4.3|l and let Ci > 0, /i G (0, 1/2). There exists Cq > such 
that if 

^ {5LB^/'^Y < Co, (A.2) 



eB5L 

then, for all z G M^, 



4 



Proof. 

Using the inequality 

1 



n>i/An> > -(pi + s)n>, 

and the result of Lemma lA.31 with q^^ = 1 — ^1, we find 

We implement unitarily the translation by z and scaling by B^^^ and end up having 
to prove that 

Z 11 
P|-1A + 1 - '^^^^ ^B{5LY-^ \B-^/^x\>^ - 2- 

But, under condition IIA.2|) this inequality is clearly true for cq sufficiently small. □ 
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